function y = uc( x)
global climit
%this returns current utiltiy of consumption level x
%note that all levels are in thousands
%so the coefficient is 1/20 = 0.05
y=1-exp(-0.05*(x-climit));
end

